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Abstract.  Two  Himalayan  ice  cores  display  a  factor-two 
decreasing  trend  of  air  content  over  the  past  two  millennia, 
in  contrast  to  the  relatively  stable  values  in  Greenland  and 
Antarctica  ice  cores  over  the  same  period.  Because  the  air 
content  can  be  related  with  the  relative  frequency  and  inten¬ 
sity  of  melt  phenomena,  its  variations  along  the  Himalayan 
ice  cores  provide  an  indication  of  summer  temperature  trend. 
Our  reconstruction  point  toward  an  unprecedented  warming 
trend  in  the  20th  century  but  does  not  depict  the  usual  trends 
associated  with  “Medieval  Warm  Period”  (MWP),  or  “Little 
Ice  Age”  (LIA). 


1  Introduction 

Relatively  little  is  known  concerning  climatic  change  over 
the  Himalayas  and  Tibetan  Plateau  for  time  scales  ranging 
from  centuries  to  thousands  of  years.  This  is  largely  due  to 
the  sparsity  of  meteorological  stations  and  its  unique  topo¬ 
graphic  setting.  The  longest  instrumental  records  from  the 
region  date  back  only  to  AD1935.  Ice  cores  collected  from 
high  elevation  glaciers  in  the  Himalayas  provided  extended 
records  of  past  climate  in  the  region  (Thompson  et  al.,  2000; 
Kang  et  al.,  2002;  Qin  et  al.,  2002). 

Ice  water  stable  isotopes  (<5lsO  and  <5D)  are  the  usual 
proxy  for  air  temperature  but  many  other  factors  may  affect 
them  such  as  atmospheric  circulation,  moisture  transport,  and 
variability  of  precipitation  (Hoffmann  and  Heimann,  1997; 
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Araguas-Araguas  et  al.,  1998).  For  instance,  the  millennial 
<5lsO  record  of  the  Dasuopu  ice  cores  suggests  a  significant 
20th  century  warming  (Thompson  et  al.,  2000)  and  existence 
of  LIA  (Yao  et  al.,  2002).  However,  paleolimnological  stud¬ 
ies  in  the  Khumbu  valley  (Fig.  1)  indicate  a  cold  phase  during 
MWP  (Lami  et  al.,  1998),  and  no  20th  century  <5lsO  enrich¬ 
ment  was  found  from  an  80  m  ice  core  recovered  in  1998 
from  the  East  Rongbuk  (ER)  Glacier  (28°0LN,  86°58'E, 
6450  m  above  sea  level.  Fig.  1)  (Qin  et  al.,  2002).  Instead, 
the  ER  ice  core  5180  record  is  significantly  and  inversely 
correlated  to  the  Southwest  Indian  Monsoon  intensity  (Qin  et 
al.,  2002;  Zhang  et  al.,  2005),  suggesting  the  dominance  of  a 
“precipitation  amount  effect”  on  stable  isotopes  (Hoffmann 
and  Heimann,  1997;  Araguas-Araguas  et  al.,  1998).  Here  we 
present  a  novel  attempt  to  reconstruct  the  summer  tempera¬ 
ture  trend  using  a  physical  proxy,  the  ice  core  gas  content, 
that  reflects  relative  frequency  and  intensity  of  melt  phenom¬ 
ena,  thus  providing  valuable  data  on  past  summer  climate 
(Herron  et  al.,  1981). 


2  Experimental  methods 

We  recovered  one  ice  core  to  bedrock  in  2001  (117.06  m 
deep,  hereafter  “Core  2001”),  and  two  cores  in  2002 
(108.83  m  deep,  hereafter  “Core  2002”,  and  95.80  m  deep, 
respectively)  on  the  col  of  the  ER  Glacier  at  6518m  above  sea 
level  (Fig.  1)  using  an  electromechanical  drill  in  a  dry  hole. 
A  repeating  survey  with  a  Sokkia  GSS1A  Global  Position¬ 
ing  System  (accuracies  of  ±5  mm  over  distance  up  to  10  km) 
in  1998  and  2002  did  not  identify  horizontal  movement  at 
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Fig.  1.  Location  map  of  ice  core  drilling  sites,  together  with  a  sketch  of  thermal  zones  on  the  glacier. 


the  drill  site.  Ice  layers  are  horizontal  in  the  cores,  suggest¬ 
ing  negligible  ice  dynamical  deformation.  The  average  an¬ 
nual  net  accumulation  is  ~400mm  water  equivalent  as  de¬ 
termined  by  snow  pits  and  the  80  m  ice  core  studies  (Hou 
et  al.,  2002).  Bore-hole  temperatures  in  the  108.83  m  core 
ranged  from  a  minimum  of  —9.6°  at  20  m  to  —8.9°  at  the 
bottom.  We  maintained  the  temperature  of  the  cores  below 
—5°  from  the  time  of  drilling  until  analysis. 


Seventy-seven  samples  from  Core  2001,  and  123  from 
Core  2002  were  measured  for  air  content  after  extracting 
the  air  by  melting  and  refreezing  the  ice  under  vacuum,  and 
chromatographic  measurement  of  the  air  peak  with  a  preci¬ 
sion  within  ±5%  (method  a).  All  samples  were  selected  be¬ 
low  the  close-off  depth  (28.2  m  for  Core  2001  and  26.2  m 
for  Core  2002,  respectively)  down  to  the  bottom  part  of  the 
cores.  Each  sample  (~50g)  corresponds  to  4-10  cm  ice  core 
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intervals  covering  seasonal  for  the  upper  core  sections  to 
multi-annual  for  the  bottom  sections.  Additionally,  we  mea¬ 
sured  4  samples  from  Core  2001  by  an  accurate  air  pressure 
method  in  a  calibrated  volume  that  has  a  proven  accuracy 
within  ±0.6%  (method  b)  (Lipenkov  et  al.,  1995).  The  ER 
results  with  method  b  are  in  good  agreement  with  chromato¬ 
graphic  measurements  (Fig.  2).  For  consistency,  we  only  in¬ 
clude  the  measurements  by  method  a  for  further  discussion. 

Gas  content  measurement  implies  that  a  portion  of  air  is 
lost  by  cut-bubbles  at  the  sample  surface.  As  the  shallower 
samples  include  bigger  bubbles,  the  cut-bubble  effect  could 
generate  an  underestimate  of  the  air  content  in  shallower 
samples  compared  with  deeper  samples.  To  account  for  this 
effect,  we  measured  the  average  bubble  size  of  six  samples 
taken  at  different  depths  between  26.4  and  115  m.  We  then 
calculated  the  cut-bubble  correction  coefficient,  ranging  from 
2%  in  the  deepest  section  to  14%  at  26  m.  We  corrected  all 
measured  air  contents  as  a  function  of  depth  with  a  third- 
order  polynomial  regression  running  through  the  six  mea¬ 
sured  correction  coefficients  (R2=0.9989). 

3  Dating 

We  dated  the  ice  core  gas  age  by  a  combination  of  methane 
and  <518Oatm  records  with  comparison  to  the  GRIP  chronol¬ 
ogy  (Hou  et  al.,  2004).  We  estimate  a  gas-ice  age  difference 
of  around  30  years  based  on  the  ice  core  bore  hole  gas  mea¬ 
surements  at  the  Dasuopu  glacier,  Himalaya  (Xu  and  Yao, 
2001),  and  calculation  from  the  present  accumulation  rate 
and  temperature  (Schwander  et  al.,  1997).  The  final  ice  age 
is  consistent  with  the  age-depth  function  of  an  80m  annually 
dated  ice  core  recovered  from  the  same  glacier  (Kang  et  al., 
2002;  Qin  et  al.,  2002).  We  estimate  a  dating  error  ±10  years 
throughout  the  whole  ice  cores. 

4  Factors  controlling  gas  content 

In  absence  of  melting,  the  air  content  of  the  ice  (V ,  in  cm3/g 
of  ice)  depends  on  the  pore  volume  ( Vc),  the  air  pressure  ( Pc) 
and  temperature  ( Tc )  prevailing  at  the  close-off  depth  interval 
(Martinerie  et  al.,  1992)  as: 


where  Pa  and  Ta  are  the  standard  pressure  (1013  mb)  and 
temperature  (273°K),  respectively;  V,  is  the  pore  specific 
volume  of  ice  at  close-off  (cm3/g). 

With  the  present  conditions  at  our  drill  site  (Pc=462mb; 
7)— 264  K),  and  the  empirical  dependence  of  Vc  on  tempera¬ 
ture  (Martinerie  et  al.,  1992),  V  is  calculated  as  0.066  cm3 /g 
using  Eq.  (1).  This  is  roughly  twice  the  average  measured  V 
values  (Fig.  2).  Thus  other  processes  during  snow  transfor¬ 
mation  contribute  to  the  low  gas  content  of  our  cores  rather 
than  those  identified  for  the  dry  snow  zone  of  polar  regions 


(Martinerie  et  al.,  1992;  Delmotte  et  al.,  1999,  and  references 
therein). 

On  our  site,  an  ice  core  annual  layer  (Fig.  1)  consists  of 
alternating  summer  percolation  and  winter  dry  snow  strata. 
In  summer,  snow  melting  can  percolate  downwards  (0° 
isotherm  at  the  end  of  summer;  Fig.  1)  (Paterson,  1994).  This 
accelerates  and  strengthens  packing  by  lubricating  the  grains. 
Refreezing  of  melt-water  also  speeds  up  the  later  stages  of 
snow  transformation.  As  a  portion  of  air  spaces  is  filled  with 
water  (Paterson,  1994),  Vc  in  Eq.  (1)  is  reduced,  resulting  in 
low  gas  content  in  the  summer  stratum. 

Energy  and  mass  exchange  experiments  carried  out  on  Hi¬ 
malayan  and  Andean  glaciers  suggest  that  snow/ice  melt¬ 
ing  is  essentially  driven  by  the  short-wave  radiation  balance 
(Aizen  et  al.,  2002;  Wagnon  et  al.,  1999).  During  the  dry 
season,  the  negative  latent  heat  flux  generates  a  strong  subli¬ 
mation  which  maintains  the  glacier  surface  cold.  During  the 
wet  season,  the  reduced  short-wave  radiation  due  to  higher 
cloudiness  does  not  compensate  for  the  increased  long-wave 
radiation.  Increased  atmospheric  humidity  leads  to  a  reduced 
humidity  gradient  between  the  snow  surface  and  the  atmo¬ 
sphere.  Thus  the  incoming  energy  generates  melting  instead 
of  sublimation. 

Because  on  our  site  the  period  with  highest  humidity  coin¬ 
cides  with  the  highest  temperature,  both  phenomena  add  on 
to  generate  melting.  Observations  on  the  Xixibangma  glacier 
points  toward  a  very  good  correlation  between  snow/ice  melt 
and  air  temperature,  and  radiation  due  to  melting  accounts 
for  up  to  86%  of  incoming  energy  (Aizen  et  al.,  2002).  Thus 
we  interpret  the  gas  content  as  primarily  an  index  of  summer 
warmth. 

We  measured  a  section  of  Core  2002  in  the  range  34.42  m 
to  35.27  m,  where  the  summer  and  winter  strata  can  be  dis¬ 
tinguished  by  seasonality  of  major  ions  (Qin  et  al.,  2002; 
Kang  et  al.,  2002,  and  unpublished  data;  Fig.  2).  The  com¬ 
paratively  low  gas  content  of  the  summer  strata  validates  the 
connection  of  gas  content  with  summer  temperature.  Winter 
layers  also  show  depleted  air  content  because  they  are  proba¬ 
bly  affected  by  water  percolation  during  the  next  summer  as 
well. 

From  Eq.  (1),  the  air  content  is  also  affected  by  changes 
of  air  temperature  and  atmospheric  pressure  (i.e.,  elevation). 
Keeping  other  factors  unchanged  under  present  conditions,  a 
±2°  change  of  Tc  results  in  a  change  of  V  of  ±0.0005  cm3/g, 
i.e.,  only  ~2%  of  the  observed  V  variations  along  our  pro¬ 
files.  Persistent  retreat  of  the  glacier  terminal  moraines  is 
evident  for  both  slopes  of  Mt.  Everest  throughout  the  late 
Holocene  (Richards  et  al.,  2000;  Zheng  et  al.,  2002).  Con¬ 
sequently,  we  expect  a  progressive  thinning  at  the  drill  site 
(i.e.,  increasing  air  pressure),  which  would  increase,  rather 
than  decrease  air  content  over  time  as  we  observe.  In  con¬ 
clusion,  we  expect  that  our  air  content  record  is  mainly  con¬ 
nected  with  the  magnitude  and  frequency  of  snow  melting  on 
the  glacier  surface. 
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Fig.  2.  Air  content  profiles  of  Cores  2001  and  2002,  together  with  SO(J  and  Na+  profiles  identifying  summer  and  winter  strata  in  the 
detailed  section.  The  sample  “Core  2002  excluded”  may  represent  an  extreme  low  temperature  period  without  any  surface  snow  melting,  and 
we  exclude  it  for  further  consideration. 


5  Summer  temperature  trend 

We  combined  the  corrected  gas  content  records  of  Cores 
2001  and  2002  into  one  time  sequence  due  to  their  simi¬ 
lar  features  (Fig.  2),  thus  leading  to  a  better  time  resolu¬ 
tion.  Both  cores  reveal  a  large  decrease  of  air  content  over 
the  last  two  millennia  from  mean  values  ~0.050cm3/g  to 
~0.028  cm3/g.  Most  of  the  decrease  takes  place  over  the 
last  200-300  years.  Air  content  decreases  in  parallel  with 


a  larger  deviation  around  the  mean  value,  which  may  reflect 
the  thinning  of  annual  layers  with  depth,  as  shallower  sam¬ 
ples  partly  depicting  the  seasonal  variability  of  air  content 
(Fig.  3a)  whereas  the  deepest  samples  integrating  this  vari¬ 
ability  through  multi-annual  depth  resolution. 

The  gas  content  profile  (Fig.  3b)  thus  suggests  unprece¬ 
dented  warm  summer  temperatures  at  the  altitude  of  our 
drill  site  for  the  20th  century.  A  comparable  warming  trend 
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Year  (A.D.) 

Fig.  3.  (a)  Average  age  representativity  of  each  sample,  (b)  Gas  content  time  series,  with  a  binomial  smoothing  trend  (thick  black  line). 
Note  that  the  gas  content  is  plotted  as  a  function  of  the  age  of  ice,  because  in  this  case  it  is  affected  by  surface  processes,  (c)  The  monthly 
resolved  Palmyra  coral  <5180  records.  The  black  horizon  represents  the  average  of  all  the  coral  dataset,  (d)  Temperature  reconstruction  of 
North  Hemisphere,  (e)  Reconstruction  of  solar  irradiance  anomalies,  (f)  Radiative  forcing  associated  with  volcanic  eruptions  recorded  in 
ice  cores,  with  a  binomial  smoothing  trend  (thick  black  line).  The  approximate  timing  and  duration  of  the  “Little  Ice  Age”  (LIA)  and  the 
“Medieval  Warm  Period”  (MWP)  are  marked  by  horizontal  bars. 


was  also  identified  in  the  Himalayas  using  meteorological 
data  (Liu  and  Chen,  2000),  dendroclimatic  signals  from  the 
Nepalese  Himalayas  (Cook  et  al.,  2003)  and  western  central 
Asia  (Esper  et  al.,  2002a).  However,  there  is  no  clear  evi¬ 
dence  in  our  record  of  other  climatic  events  such  as  the  “Lit¬ 
tle  Ice  Age”  (LIA)  or  the  “Medieval  Warm  Period”  (MWP), 
thus  it  does  not  reinforce  the  trends  based  on  the  reconstruc¬ 
tion  of  Northern  Hemisphere  temperature  for  the  past  two 
millennia  (Fig.  3d;  Mann  and  Jones,  2003),  or  the  dendro¬ 
climatic  signal  from  the  Northern  Hemisphere  extratropics 
(Esper  et  al.,  2002b).  On  the  other  hand,  a  palaeoclimate  re¬ 
construction  based  on  pigment  and  diatom  assemblage  varia¬ 
tions  from  two  lakes  of  the  Khumbu  valley  (Fig.  1)  indicates 
warm  phase  during  the  13— 16th  centuries  and  cold  phase  dur¬ 
ing  the  5-12th  centuries  (Lami  et  al.,  1998).  The  relatively 
lower  temperature  before  the  12th  century  is  consistent  with 


the  general  advance  of  Himalayan  glaciers  during  the  5-6th 
and  9-12th  centuries  (Rothlisberger  and  Geyh,  1985).  The 
recent  coral  reconstruction  from  Palmyra  Island  in  the  cen¬ 
tral  tropical  Pacific  also  supports  a  cool  MWP  and  a  warm 
LIA  (Fig.  3c;  Cobb  et  al.,  2003),  which  indeed  better  fit  with 
the  trend  from  our  gas  content  record. 

Though  decadal  to  centennial  climate  variability  during 
the  pre-industrial  appears  to  have  been  dominated  by  solar 
and,  to  a  less  extent,  volcanic  forcing  (Shindell  et  al.,  2003), 
connection  between  our  signal  and  solar  (Fig.  3e;  Bard  et  al., 
2000)  or  volcanic  (Fig.  3f;  Crowley,  2000)  activity  variations 
for  the  past  millennium  is  unclear.  The  primary  impact  of 
solar  activity  on  climate  is  to  change  thermal  gradients  rather 
than  direct  solar  heating.  This  means  that  the  high  and  per¬ 
haps  mid-latitudes  could  be  out-of-phase  in  temperature  with 
the  low  latitudes  because  of  displacement  of  heat  via  changes 
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in  atmospheric  and  oceanic  thermohaline  circulation  (Goosse 
et  al.,  2004).  The  resemblance  of  our  summer  temperature  re¬ 
construction  (Fig.  3b)  and  the  tropical  Pacific  corals  (Fig.  3c) 
during  MWP  and  LIA  might  reflect  interaction  between  the 
El  Nino  -Southern  Oscillation  (ENSO)  and  the  Asian  Mon¬ 
soon  (Charles  et  al.,  1997).  The  mechanisms  causing  such 
variability  may  involve  the  large-scale  atmospheric  circula¬ 
tion  encompassing  the  Pacific  and  Indian  Ocean  monsoon 
areas.  Our  climatic  signal  suggests  that  such  a  spatial  pat¬ 
tern  might  exist  before  the  involvement  of  modern  human 
activities  in  the  tropical  climate  system. 

6  Conclusions 

Though  a  quantitative  calibration  to  temperature  is  unavail¬ 
able  due  to  lack  of  ad  hoc  meteorological  data,  coarse  sam¬ 
pling  resolution  and  uncertainty  of  ice  core  dating,  our  gas 
content  reconstructions  highlight  an  unprecedented  warming 
trend  in  the  20th  century  compared  with  the  last  two  millen¬ 
nia.  Moreover,  it  suggests  that  Himalayan  climate  might  be 
out  of  phase  with  extratropical  climate  on  centennial  scale, 
highlighting  the  complexity  of  climatic  change  in  the  tropic. 
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